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ABSTRACT 

We generahze previous calculations to a fully relativistic treatment of 
adiabatic oscillations which are trapped in the inner regions of accretion 
disks by non-Newtonian gravitational effects of a black hole. We employ the 
Kerr geometry within the scalar potential formalism of Ipser and Lindblom, 
neglecting the gravitational field of the disk. This approach treats perturbations 
of arbitrary stationary, axisymmetric, perfect fluid models. It is applied 
here to thin accretion disks. Approximate analytic eigenfunctions and 
eigenfrequencies are obtained for the most robust and observable class of modes, 
which corresponds roughly to the gravity (internal) oscillations of stars. The 
dependence of the oscillation frequencies on the mass and angular momentum 
of the black hole is exhibited. These trapped modes do not exist in Newtonian 
gravity, and thus provide a signature and probe of the strong-field structure of 
black holes. Our predictions are relevant to observations which could detect 
modulation of the X-ray luminosity from stellar mass black holes in our galaxy 
and the UV and optical luminosity from supermassive black holes in active 
galactic nuclei. 
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1. INTRODUCTION 

Our group has been investigating consequences of the reahzation [originally by Kato & 
Fukue (1980)] that general relativity (or similar relativistic theories of gravitation) can trap 
normal modes of oscillation in the regions of accretion disks near black holes (or weakly 
magnetized neutron stars). Since no such trapped modes exist in Newtonian gravity, they 
provide a signature for the existence of strong gravitational fields. The frequencies of these 
modes serve as probes of the structure of the field as well as the physical conditions within 
the accretion disk. 

Using a modified Newtonian potential, adiabatic eigenfunctions and eigenvalues of the 
lowest acoustic (p) modes (Nowak & Wagoner 1991, hereafter NWl) and internal gravity 
(g) modes (Nowak & Wagoner 1992, hereafter NW2) have been calculated. Their rates 
of growth or damping due to anisotropic (turbulent) viscosity and gravitational radiation 
reaction were also determined (NW2). Finally, the amplitude and coherence time of the 
oscillations, and the resulting fractional modulation of the spectral luminosity of the disk, 
have been estimated (Nowak & Wagoner 1993, hereafter NWS). Although this modulation 
appears to be at most a few percent, we feel that the potential payoff certainly justifies 
dedicated monitoring of a variety of supermassive black hole (AGN) candidates as well as 
the appropriate binary X-ray sources in our galaxy. 

In this paper we begin to reanalyze such modes of oscillation in terms of perturbations 
of the general relativistic equations of motion of perfect fluids within the Kerr metric. Ipser 
& Lindblom (1992, hereafter IL) have developed an elegant formalism to analyze small 
adiabatic pulsations of relativistic perfect fluids which in equilibrium are in pure rotation 
in stationary and axisymmetric spacetimes. We apply this formalism to a thin accretion 
disk around a Kerr black hole. In this case, the Cowling approximation of negligible 
self-gravitational effects of the disk should be valid. All perturbed quantities can then be 
obtained from a scalar potential governed by a single linear two-dimensional second-order 
partial differential equation. With the appropriate boundary conditions, this equation 
yields trapped modes with discrete spectra which in combination can determine the black 
hole mass M and dimensionless angular momentum parameter a = cJ/GM'^. 

We then focus on the analogs of the gravity modes that exist within stars. They should 
be the most observable trapped modes since they appear to occupy the largest area of the 
disk, including its maximum temperature. Approximate expressions for their eigenfunctions 
and eigenfrequencies are obtained, and explicit results are presented for the case of small 
buoyancy frequency. 
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2. PERTURBATIONS OF BLACK HOLE ACCRETION DISKS 

2.1. General Relations Within Thin Disks 

Unless otherwise indicated, we will adopt the conventions and notation of IL (such as 
gravitational units: G — c — 1). In addition, all lengths and times are made dimensionless 
by expressing them in units of the black hole mass M. The matter within the accretion 
disk is assumed to be represented by a perfect-fluid stress-energy tensor: 

= (2-1) 

where p is the mass-energy density, p is the pressure, is the four-velocity of the fluid 
(with M^Ma = — 1), and the operator q"-^ = g""^ + u°'v}' projects components perpendicular 
to u"-. In the neighborhood of the equatorial plane 2; = 0, we use the Novikov and Thorne 
(1973) expression for the Kerr metric in Boyer-Lindquist coordinates : 

r^A* A* 

ds^ = —dt^ + ^id(j)-ujdt? + —dr'^ + dz^ (2-2) 

^^^^^ 

(dV _ _j^(d_ r^(d_y ^(d_y (d_y 

[dsj ~ r^A^yat^'^dcp) ^A*[d(p) [dr J ^ [dz J ' 

with 

A* = r^-2r + a% A* = + r V + 2ra% u; = 2ar/A\ (2-4) 

where corrections of order [z/rY are neglected (but not later, in the vertical equations of 
motion). Note that \\gab\\^ -r'^ and gu = -A*g't"t' = -1 2/r. 

We further assume that \u^\ <^ \u^\ <C \ru'^\ within the stationary and axisymmetric 
unperturbed disk, in which the properties of the fluid are also symmetric about the 
equatorial plane. The fluid particles are thus travelling in approximately circular orbits 
(assumed in corotation with the black hole) with four- velocity 

= /3(t» + 0(^») . (2-5) 

With x"- — {t,r,(t),z}, the Killing vector components are simply t"' — (1,0,0,0) and 

= (0, 0, 1, 0). Within the unperturbed disk the Euler equation, in the form Qac^bT^'^ — 0, 
becomes 

u'^bUa = -{^aP)/P + P<p\Van = -(V„p)/p , (2-6) 

where Va is the covariant derivative. We have set p -|- p ~ p since the speed of sound 
Cg ^ hQ <^ rQ < c — 1, where h{r) is the effective half-thickness of the disk. 
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The functions (5 (equivalent to the function 7 employed by IL) and Q (the angular 
velocity measured by an observer far from the disk) are well approximated by the 
free-particle expressions 

3^ + » n= ^ (2-7) 

P ^3/4(^3/2 _ 3^1/2 + 2a)V2 ' " r3/2 + a' ^ ^ 

since the disk is assumed thin. Away from the inner edge of thin, and even slim 
(Abramowicz et al. 1988) disks, the effects of the actual radial velocity and radial pressure 
gradient on our eigenmodes will be small (NW3). 

Next, we consider the linearized equations for the Eulerian (at a fixed point in 
spacetime) perturbations 5p, 5p{<^ Sp), and 5u"'. We obtain these by perturbing the two 
projections of the energy-momentum conservation equation, the law of energy conservation 

UbVaT''''^0 O^U^VaSp + Su^VaP + pVaSu", (2-8) 

and the Euler equation of motion 

qy^j^bc^Q ^ = p(^q^u^^M'' + Su^^bu'') + 6pu''VbU'' + q''''VbSp, (2-9) 

where Su"- — q^Su^ and UaSu"- — 0. 

We impose the adiabatic condition that the perturbations preserve the entropy per 
baryon (s). This condition implies the following relationship between Lagrangian variations 
(associated with a given element of fluid) in the pressure and the density: 

Ap = (rp/p)Ap . (2-10) 

Although IL generalized V to be an arbitrary stationary and axisymmetric function in the 
equilibrium fluid, for these adiabatic perturbations V = {p/p){dp/dp)s = {p/p)c^g. Under 
most conditions (gas or radiation pressure dominance) , the equation of state will be of the 
form p — K{s)p^, so we shall take the adiabatic index F (denoted 7 by NWl-3) to be at 
most a slowly varying function of r. 

One can relate the Lagrangian and Eulerian variations via the Lagrangian displacement 
vector giving Ap = Sp + i°'VaP and Ap = 5p + ^"'VaP- The negative Lie derivative 
of the unperturbed velocity gives its Eulerian perturbation (Friedman & Schutz 1978): 
5u°' — Qbiu'^'^ci^ — ^"^VcU^). The previous equation imphes a gauge freedom in 
corresponding to the addition to ^" of any function parallel to u"" (Schutz & Sorkin 1977). 
We choose ^°-Ua = to fix the gauge. 

From now on, we shall work with the azimuthal and temporal Fourier components 
of all quantities such as the displacement vector ^" — > ^"(r, z) exp[i(m0 -\- at)]. We also 
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employ the corotating frequency uj{r) = a + mQ{r). (Note that IL exchange the definitions 
of uj and a; we employ those of NWl-3.) This Fourier decomposition provides us with an 
algebraic relationship between and Su"". 

r = -tiuj(3)-^[q^, - (z/a;)g>^Vbfi](5^i'' . (2-11) 



With p <^ p as usual, the scalar potential employed by IL becomes 6U = 6p/ p, the 
specific adiabatic Eulerian perturbation of the enthalpy. However, we will also find it 
convenient to use their rescaled potential 

6V = 6U/{l3uj) = 6p/{p(3uj) . (2-12) 

Following IL, one can then obtain the relation 

(2-13) 



in terms of the function 



pT ujj3 



Aa = ^-^ = p-'VaA , (2-14) 

ppl 

where the final equality is valid if F is constant, with A = \n{p/p^/^)+ constant. Note [from 
equation ( |2-10|) ] that Aa vanishes for any barotropic equation of state p = p{p). 

Since the accretion disks that we shall consider are thin, vertical gradients of 
equilibrium quantities (such as pressure) are typically much greater than their radial 
gradients. Therefore the Brunt- Vaisala or buoyancy frequency, A^^, is given by 

A^VaP = A'V.p = P'^Nl . (2-15) 



Radial perturbations of the circular orbits of free particles (corresponding to a 
zero-pressure, infinitesimally thin disk) oscillate at the epicyclic frequency k given by 
(Okazaki, Kato, & Fukue 1987) 

2u^n, - 2u^n, = p'k' = p'n' (i-^ + l^^^) , (2-16) 

where u"" (not to be confused with the corotating frequency uj) and fi" are the vorticity and 
angular velocity vectors (IL). Similarly, vertical perturbations of free particle circular orbits 
produce oscillations at the vertical frequency f2_L given by (Kato 1990) 
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The frequency Q± also appears in the equation of vertical hydrostatic equilibrium (Kato 
1993): 

dp/dz = -pfS'^nlz , (2-18) 
and therefore in our definition of disk half-thickness h(r): 

pir,0)/p{r,0) = {hf3nj_y . (2-19) 

We note that equation ( |2-18| ) is incorrect in Novikov & Thorne (1973), as also pointed out 
by Riffert & Herold (1995). We can now obtain the useful expression 

= -Ql{zdA/dz) , (2-20) 

from equations ( p-14|) and ( |2-15| ). 

These frequencies (k, Q±, and A^^) play the key role in determining the modes of 
oscillation of the disk. In Figure |1|, we show the radial dependence of the free-particle 
orbital frequencies Q, k, and Q±, for three values of the angular momentum parameter 
a. Thorne (1974) has shown that a black hole can be spun up by accretion to a value of 
a = 0.998, only slightly less than its absolute limit of unity. 

The equations which govern the perturbations of the disk simplify considerably if the 
perturbations of the gravitational field are negligible (Cowling approximation), as is true 
for thin accretion disks. Using the perturbed equation of motion ( |2-9| ) , IL are then able 
to algebraically relate the Eulerian perturbations 6u"', 6p, and 6p to 6U and its covariant 
derivative. Using the perturbed equation of energy conservation (p-8|), IL then obtain a 
two-dimensional partial differential 'master' equation [their equation (39)] which governs 
the rescaled potential 6V. For our problem it assumes the form 



{A*y/^pH^"'DJv] + {A*y/^uj(3^SV = , (2-21) 



where Z)^ {p = r, z) is the covariant derivative associated with the metric g^'^ of the plane 
orthogonal to the Killing vectors. We find (from the expressions in IL) that 

= g--uj\uj^ - , H'^ = uo\uo^ - Nl)-\ (2-22) 

and W'^ is negligible. In addition, we obtain 

$ = puPcf - p\ujp)-^H''Al + W.[p\ujp)-^H''A,] . (2-23) 

We have restricted consideration to low azimuthal modes: |m| <^ r/h. The reduction of 
the function $ to the above form also depends upon the assumption that the corotating 
eigenfrequency uj lies in the range h/r ^ Ic^l/fi ^ r/h. 
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2.2. Radial- WKB Approximation and Separation of Variables 



IL applied their potential formalism to two problems: general adiabatic perturbations 
of nonrotating stars and barotropic perturbations of rigidly rotating stars. As in NW2,3, we 
shall apply it to thin disk perturbations whose radial wavelength A(r, z) = 2'k\6V {d6V / dr)~'^\ 
is much less than r, the scale of radial variation of equilibrium properties of the disk. This 
turns out to be a good approximation, and it was also found in NWl-3 that A ~ h{r) for 
most of the lowest modes studied. Therefore, we work to lowest nontrivial order in the 
WKB parameter A/r. The master equation (p-21| ) then becomes 



d 
or 



d5V 



dr 



+ - — 

p dz 



P 



dSV 



{uj^-N]) dz 



+ 



1 d 



P 



c2 p^^dz {{uj^-N^) dz 



6V = 0. 



(2-24) 



We have also allowed for rapid variations near the radial corotation resonance at ui"^ 



K 



In the quasi-Newtonian limit, NW2 neglected the terms involving dlu"^ 
corresponding equation (2.4). 



Nl)/dz in their 



In order to express Lagrangian perturbations such as the boundary condition 
= ujl3p5V + aP = 0, we need the expressions 



d 



dr 



6V 



(3(00^ - A2) 



dz 



5V + pAJV 



for the relevant components of the displacement vector. We also find that 



where u% = -f3^[l - (2/r)(l - a£l)\. 



-u-ut) 



ruj 



'dn 



(2-25) 



(2-26) 



As in previous work (NW2,3), equation ( P-24| ) is separable in r and rj = z/h{r) to 
lowest WKB order, if we assume that c1{r,z) = c'^{r,0)g{r,ri) is a slowly-varying function 
of r. We note, from our definition ( p-19| ) of h{r), that c^(r, 0) = r(/i/5f2_L)^. In order to 
bring the radial separated equation to a familiar form [cf. short wavelength perturbations 
of nonrotating stellar atmospheres (Hines I960)], we take the slowly varying separation 
function to be 



T{huj) 



"^(r) - 



The redefined separation function \E'(r) was denoted rT(r) by NW2,3. So now we set 
6V = Vr{r)Vri{r,ri), with Vrj a slowly-varying function of r, and obtain the separated 
equations 



d^Vr 

dr'^ 



^ _ ^2)^ ^ ^2(^2 _ ^2)(^ _ ^/^2)Vr = , (2-27) 



dr 



- 8 - 



where tu* = u /Vt_[_ and 



dA 

dr] 



7] dA^ 



+ 



c2(r,0) Thmi ' 



^* df] J dr] 



dr] 

d^A ^dA 
dr] 



K = 



7^ = 



fi2 



(2-28) 



(2-29) 



Equations (|^ and (|^) agree with equations (2.3a) and (2.3b) of NW3 in the 
quasi- Newtonian hmit, except for the term proportional to dluj"^ — n'^)/dr missing in NW3. 



2.3. Boundary Conditions 

We shall take the top and bottom surfaces of the disk [r] = ±r]o{r)] to be defined by 
the vanishing of the equilibrium pressure p. The boundary condition on the perturbations 
there is taken to be the usual one: = Ap = Sp + zP- Using equations ( P-12| ), ( P-18| ), 
( F20| ), (p^ , and ( p^ ), this becomes 

where B{5V) = lo^SV - r]d5V/dr]. 

We shall also make the physically plausible assumption that the buoyancy frequency 
A^^ is finite [vanishing on the midplane according to equation ( |2-20D ], and remains so as 



p ^ at the boundary. It then follows that dA/dr], and therefore also A and d'^A/dr]'^, are 
finite. From hydrostatic equilibrium, we then obtain limiting expressions for the following 
unperturbed functions: 



p-f{r){Vo-v)^^, 9=^^)--{^)voiv-Vo) (2-31) 



as ?7 ±?7o 



With this behavior near the boundary, the vertical separated equation ( p-28| ) assumes 



the form d'^Vrj/dr]'^ + B{yr,)/{Tg) — >• 0. However, the boundary condition (p-30|) is satisfied 
if pBiVri) vanishes, while the stronger condition that Ap/p be finite requires that B(yrf)/g 
remains finite at the boundary. These conditions thus translate into the basic requirement 
that d'^Vrj/dr]'^ be finite at the boundary. Since the boundary is a singular point of the 
first kind for equation ( P-28| ) (see also section 3.3), it can be shown that one of the two 
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independent solutions of this equation is regular at r] 
finite there. 



ir^o, with equivalently being 



It is then seen from equations (|2-28|) and (p-29|) that in order that these regular 
solutions remain finite through the disk, we must require that remain less than u;^ to 
avoid the 'buoyancy resonance' at uj'^ = N'^ which occurs in equations (|2-24|) and (|2-29|). 



3. G-MODE ANALYSIS: EIGENFUNCTIONS 

We are now in a position to attempt to obtain analytical estimates of properties of 
the relativistically trapped modes of oscillation of the disk. We will concentrate on the 
analog of the gravity (g)-modes, since they appear to be least sensitive to the poorly known 
physical conditions at the inner (and outer) radial boundary of the disk. In addition, 
they typically occupy the largest area of the disk. The pressure (p)-modes and (m = 1) 
corrugation modes usually occupy a narrow region bounded by the inner edge of the disk 
(NW2,3; Perez 1993; Kato 1990; Ipser 1994). 

The g-modes are found to exist for values of the separation function \E' which satisfy 

^>ujI [uj, = uj/n^], (3-1) 



corresponding to /t^ — cu^ > from equation ( 2-27|) . Relativistic gravity forces the radial 



epicyclic frequency k to reach a maximum kq at a radius tq which is larger than that 
of the inner edge of the disk r = r, (very near the innermost stable circular orbit). 
Therefore we can require that these modes become effectively evanescent at r < T_(cr) and 
r > r^{a), where r_ and r_|_ (> r_) are the places where fi;^(r) = a;^(r). Recall that the 
eigenfrequencies a (measured far from the disk) are related to the corotating frequencies 
by cr = uj{r) — mQ{r), and that Q{r) > K{r) > 0. In Figure ^ we plot these functions that 
determine the radii r±{a) between which the eigenfunctions are concentrated. In Figure |^ 
we show the dependence of these boundaries of the trapping region on the eigenfrequency. 

We see that for m = these g-modes exist if cu^ = cr^ < Kq. For m > they exist 
if > cr > — max{mQ{r) + ^(r)} = — (mi^^ + Km), with the maximum occuring at 
r = Vm- Since only the square of u appears in our fundamental equation (|2-24[) for the 



scalar potential 6V, the case of negative m is obtained by the replacements m — m, 
o" ^ —cr, as is also seen from Figure Therefore since only |cr| is observable we can 
assume without loss of generality that m > and a < throughout this paper. Thus 

< — cr < (Trn = K,m + TTl^lm- 
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The eigenfrequencies are also bounded by the condition 

max{N^} < Lu^ . (3-2) 

This inequahty follows from the buoyancy resonance avoidance discussed at the end of the 
previous section. The maximum refers to the behavior of N"^ as a function of z (for each 
value of r) and is typically expected to be reached at the boundary of the disk. 



3.1. Radial WKB Solutions 



Within our radial WKB approximation, the radial separated equation (|2-27|) can be 
written in the form of the Schrodinger equation 

d^W/dr^ + U{k^ - uj^)W = , (3-3) 

where W = {k^ — u^)~^dVr/dr and 11 = a'^[{'^/ul) — 1]. From equations ( p^-25| ), we see 
that the displacement components C,^ oc W and oc K = —Tlr^dW/dr. In the vicinity 
of r = r^(/i = — , +), we have YI^h? — uj"^) = Yl^V^ir — r^) if is not close to r^. (We 
will treat the case of close to later.) Here the various quantities = ^iid 
V = d{H? — uj^^jdr. [Recall that uj = uj{r) is a constant only if m = 0.] 

Near r+, we introduce the coordinate x+ = g+(r — r_|_), g_|_ = (n|P|)y^. This transforms 
equation ( p-3D in the vicinity of r_|_ into the Airy equation d'^W/dx'^ — x+W = 0. Its solution 
which exponentially decays for r > r_(_ is 

W cx Ai{x+) = n-\x+/3y^'^Ki/3{2x^f /3) , (3-4) 

where K,y{X) [and Iu{X)] are the modified Bessel functions of imaginary argument. For 
g_i_(r+ — r) ^ 1, this has the asymptotic form 

ly oc (r+ - r)^^/^ cos[(2/3)g+ ^(r+ — r)^/^ — n/A] . (3-5) 

Near r_, the governing equation is the same with x+ replaced by x_ = q-{r^ — r), 
q_ = {Il\V\)^J'^. However, because the inner edge of the disk (at r^) is not far from r_, we 
must now include the other independent Airy function 

Bz{x^) = (a;_/3)^/2[/_i/3(2/_/V3) +/i/3(2xi/V3)] , 



so that in the vicinity of r_ 

W oc Ai{x^) + eBi{x^) . (3-6) 
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The mixing coefficient 9 of tlie part tliat grows exponentially for r well inside r_ is 
determined from the boundary condition at rj. 

The asymptotic forms of this solution are 

ly oc (r - r_)-^/^ cos[(2/3)g^/^(r - r^)^/^ + (tt/A) - arctan(l/2e)] (3-7) 

for g_(r — r_) ^ 1, and 

ly oc (r_ - r)-^/^{exp[-(2/3)g'/'(r_ - rf/^] + 20 exp[+(2/3)g'/'(r_ - r)^/'']} (3-8) 

for q-{r- — r) ^ 1. Away from r_ and r+, we may employ the standard WKB 
approximation, which takes the form 

W oc - cu2)]-i/4 cos (^j'^ [II{k^ - uj^)Y'^dr' + $oj (3-9) 

well within the interval {r_,r_|_}, and 

W oc [n(c^2 _ ^2)]-i/4 ^Qgj^ (j^'[n(cu' - /s:2)]^/'cir' + $3) (3-10) 

well inside r_. The phases $fc should not be confused with the function $ given by equation 
( ^-23| ). We also define the phases 

We now match these WKB solutions to the appropriate asymptotic forms of our local 
solutions. 

The first overlap region exists at values of r slightly less than r_, but such that 
(r_ — r) ^ l/o'-- We first note that the WKB expression ( P-IQ ) assumes the form 

oc (r_ - r)-^/^cosh[(2/3)g!/^(r_ - rf/^ + $1 + $3] 



as r — i> r_ from below. Comparing this with equation ( |3-8| ) , we see that these solutions 
match if 

$1 + $3 = |ln(20) . (3-11) 

The second overlap region exists at values of r slightly greater than r_, but such that 
(r — r_) ^ l/o'-. There the WKB expression ( p-9| ) assumes the form 

ly oc (r - r_)-^/^ cos[(2/3)g!/^(r - r^f^ + $0] 
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ly oc (r+ - r)-^/^ cos[(2/3)g+/^(r+ - r)^/^ + $q + $ 



as r — » r_ from above. Comparing this with formula ( P-7|) , we see that these solutions 
match if 

$0 = (7r/4) - arctan(l/20) + nin , (3-12) 

where ni is an integer. The final overlap region exists at values of r slightly less than r+, 
but such that (r+ — r) ^ l/q'+- There the WKB expression ( |3-9| ) assumes the form 

2] 

as r — > r+ from below. Comparing this with formula ( |3-5| ) , we see that these solutions 
match if 

<l>o + $2 = -(7r/4)+^27r, (3-13) 

where n2 is another integer. 

Combining equations ( |3-11|) , ( P-12|) , and ( |3-13| ) then gives us our first fundamental 
WKB relation between the eigenfrequency a = uj{r) — mQ{r) and the separation function 
\E'(r): $2 — arctan{exp[— 2($i -|- $3)]} = {n + 1/2)tt, where n is also an integer. Although 
$1 and $2 are equal to positive definite integrals of known functions of the unperturbed 
disk, $3 is related to the still unspecified mixing coefficient 6 by relation ( p-ll| ). 

Note that on the midplane z = 0, Ap = 6p + EJdp/dr. Assuming that the first term 
dominates, requiring that Ap/p remain finite as r — > means that dW/dr = —UVr —>■ 0. 
Applying this boundary condition to equation ( P-10| ) then requires that $3 = 0, equivalent 
to 6* = |exp(2$i), and therefore 

$2 - arctan[exp(-2$i)] = {n + l/2)7r . (3-14) 

This is the first fundamental (implicit) WKB eigenfrequency relation. 



3.2. Local Approximation 

We will see that the lowest g-modes (those which have the smallest number of radial 
nodes in their eigenfunction, corresponding to the highest frequencies) are concentrated in 
the neighborhood of r = rm- The largest value of |cr| occurs where r_(cr) and r+(cr) merge 
at r = rm, which is where u = — /t and (i(K^ — uj'^)/dr = 2nm{dn/ dr + mdfl/dr) = 0, as 
shown in Figure ^ Recall that m > 0, a < 0, and quantities evaluated at r = (and, if 
necessary, a = —am) have a subscript m, with the maximum allowed value of |cr| denoted 
by am = Km + mVtm. 

Thus it is reasonable to look for an approximate analytic solution W{r) to the 
reformulated radial equation ( ^-31 ) via a local approximation of the function K^(r) — co'^(r). 
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Therefore we employ the expansion (truncated at first order in |cr| — cTm and second order 
in r - Tm) 

K^-uj'' = 2K^{a + am) - V2{r - r^f , (3-15) 

where 1^2 = ~(l/2)[c?^('t^ — '^'^)ldr^\m- We ignore the less critical radial dependence of 
n = a^i^jJl - 1) in equation (|3), letting n(r) = H^. 

If we also employ a dimensionless radial distance x and dimensionless eigenvalue A 
defined by x'^ = IlmT^2ii" — rm)^ and A = 2(nm/P2)^^^'«m(o" + o"m), equation (|3-3|) assumes 
the form of the quantum harmonic oscillator equation (PW/dx"^ + (A — x^)W = 0. Its 
eigenfunctions which decay to zero far from and the corresponding eigenvalues are given 
by the well-known expressions 

Wn oc H^{x) exp{-lx^) , A = 2n + 1 (n = 0, 1, 2, . . .) ; (3-16) 

where if„ are the Hermite polynomials. They may be constructed from the expression 
Hn = (— l)"'exp(x^)[o?"'exp(— x^)/(ix"]. We shall later show that these eigenvalues A„ are 
the same as those obtained in this local limit from the general expressions of the previous 
section, so we have used the same index n. 

Thus we obtain the first-order expression 

am - K\ = iV^/IlmY'^K-;^{n + 1/2) (3-17) 

for the splitting of the eigenfrequencies an- Note, however, that we have not yet included 
the splitting produced by the vertical eigenfunctions. To obtain a first rough estimate of 
the magnitude of the frequency splitting, define an effective radius of curvature Rm by 
1)2 = K^l and use Cg ~ hVt ~ hVL± in the expression ( |2-29D for in 11. We then find 
that 

<y^ - ~ (n + l/2)VLm{h/R)m[^nl/t^^ - 1^'^ • (3-18) 

Thus the lowest modes of thin disks indeed have eigenfrequencies very close to a^- On the 
plots in Figure |^ we have exhibited the dependence of this key frequency on the black hole 
angular momentum parameter, for m = 0, 1,2,3. (As with all orbital frequencies, am is 
inversely proportional to the mass of the black hole.) 



3.3. Vertical WKB Solutions 

We next attempt to obtain, in a way similar to that employed in Section 3.1 for the 
radial separated equation, solutions to the vertical separated equation. It has the form 



d^Vjdif + ai(r/,r)9K,/9r/ + a2{r],r)Vr, = , 



(3-19) 
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with the expressions for the coefficients given in equation ( 2-28|) . Here we shall also employ 
a WKB approximation, which in this case will be valid if 

\d\nVrf/dri\ ^ \d\na2/dri\ ~ \d\nai/dri\ ~ l/{r]o — rj) . 

This condition is not as easily satisfied as in the radial case, since the properties of the 
unperturbed disk (in ai and 02) vary much more rapidly in the vertical direction. However, 
we may take the separation function \1/ (in 02) to formally be the required large WKB 
parameter. Then the standard WKB analysis gives (away from the boundaries rj = ir^o) 
the approximate solution 



T/ -1/4 

V,j (X a2 exp 



cos 



(3-20) 



ai{r]')dri' 

where X = (or 1) for an even (or odd) function of 1]. 

In our previous analysis of the behavior of the equilibrium properties of the disk near 
its boundary, we isolated the dominant terms, which produce the following asymptotic 
behavior of the coefficients in equation ( p-19|) : 

-l/[(r - l)(r/o - v)] , «2 ~ ^.V[^o(r - l)(r/o - v)] 



ai 



as ?7 — > ±?7o. Integrating these functions then gives us the form of the WKB solution ( |3-20D 
near the boundary: 



oc (r/o - r])-W^ cos [q{r]o - i^f'^ - Bq + X7r/2 
where q = 2uj,[r,QiT - l)]-^'^ and 60 = /J"' J^^dr]. 



(3-21) 



Next, we obtain the behavior of the exact solution of equation (|3-19| ) near the boundary 
(?7 — i> rjo), using the above limiting expressions for ai{ri) and 02(77). The nonsingular solution 
to the limiting equation is then found to be 



(2-r) 



(3-22) 



involving a Bessel function of the ffist kind. Now if (r/o — 1]) ^ 1/q'^, it assumes the same 
form as the limiting WKB solution (|3-21|) , but with a specific phase shift 7r(3 — r)/[4(r — 1)]. 
Matching the two expressions gives the second fundamental (implicit) WKB eigenfrequency 
relation 

rrio I r 1 z.'^ — r\ Tl 

(3-23) 



60 



a2{r])dri = vr 



i + 7 



X 
2 



4 vr - 1 

where j is an integer. We note that with Qq positive, j > —(3 — r)/[4(r — 1)] for the even 
eigenfunctions (X = 0) and j > —(1 + r)/[4(r — 1)] for the odd eigenfunctions (X = 1). 

As usual, all the results in this section hold for any axial mode number m, which 
appears in the representation of the corotating frequency u!{r) = cr + mf2(r). 
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4. EIGENVALUES AND EIGENFREQUENCIES 



In our WKB approximation, the vertical eigenvalues \E'j(r, cj^) of the separation 
function and the corotating eigenfrequencies uj = uJnjij) = anj + mVL{r) are obtained as 
roots of the radial and vertical transcendental equations ( p-14 ) and (|3-23|) , respectively. We 
reproduce them in the following forms: 



a[{^/u:l - - uo^p^dr - arctan exp -2 / " a[{^ /ul - l){uo^ - K^)Y/Hr 



(n + l/2)7r 



- ujI uldp/df] 

r r?7p 

where we require n > and j + Sy > 0, with 

3-r 



+ Q 



1/2 



dr] = (j + 6r)iT 



5r 



4(r-l) 

We have also introduced the quantity 



(even modes) 



1 + r 
4(r-i) 



OA 

dr] 



(odd modes) 



(4-1) 
(4-2) 

(4-3) 



(4-4) 



Note that a^(r) = P'^grrPc/iXPc)- Recall the restriction (|3-1|) on the separation function \Efj 
and the lower limit ( p-2|) and the upper limits on the eigenfrequencies given in the previous 
section. 



Equation ( |4-1| ) implies that 



{n + l/2)n < / a[{^j/ujl - 1){k^ - uj^)Y/^dr < {n + 1)71 . 



(4-5) 



Since for any given value of j, \E'j(r, a;) is bounded everywhere in the (finite) interval 
r_ < r < r_|_, as are all the other quantities in the above integrand, the first inequality 
shows that n has a finite maximum value. Thus for any given value of j, there exist only a 
finite number of eigenfrequencies anj- On the other hand, equation ( [4-2| ) then shows that 
00. Therefore, for any given value of n > 0, the second inequality in 



\E'j 00 for j 



equation ( [4-5| ) implies that |cr,„j| —>■ Hm + mQm = (for any m > 0) from below (with 
r± ^ ± 0) as j 00. 

Using the fact that the left hand side of equation ( |4-1|) is a continuous function of u it 
follows that for a given value of n this equation has at least one root anj for any sufficiently 
large value of j. We can therefore summarize the above results as follows. 
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The spectrum of eigenfrequencies anj in the range — mf2(rj) > a > —[mQ{rm) + i^{rm)] = 
—dm is discrete, with the only condensation point Om- For any given integer j satisfying 
j + 5r > 0, the number of eigenfrequencies is finite, with n = n£)„, n£]„ + 1, . . . , n^)j^. For 
any given integer n > 0, there is a sequence of eigenfrequencies which satisfies 

hm anj = -{Km + rriQrn) = -o-m ■ (4-6) 

For 771 = 0, either sign is allowed, but for convenience we have chosen the negative. Typically, 
there will be monotonicity in n and j; that is, |cr„+ij| < \(Tn,j\ and > 



4.1. Local Approximation 



For eigenfrequencies close to am, with K^(r) — to'^(r) approximated by equation ( ^-15| ), 
the first term on the left hand side of equation ( [4-lD becomes small as r_,r+ — > r^, but 
the second one is exponentially small and thus may be neglected. After some algebra, this 
equation then provides the approximate expression 



^2 



(Jm. — \(y, 



Kma{rm) 



-1/2 



71 + 



(4-7) 



This expression is consistent with the previous quantization ( |3-16| ) of the eigenvalue A„ 
obtained from a more specific analysis of the local approximation, as can be seen from the 
definition of A in that section. 



Using this expression and the fact that \E'j oc for large j from equation ([4-2|) , we see 



that am — I'^njl oc (n + 1/2)/ j for j — >• oo. 



4.2. Small Eigenfrequencies (for m = 0) 

Let us now investigate the low frequency regime of the radial (m = 0) modes, in which 
r_ — i> Ti and r+ — oo, which is in a sense opposite to that of the local approximation. As 
uj = a 0, the first term in equation ( |4-1| ) becomes 

/ a{r)[{^/ujl-l){K^ -uj^)f^dr^\uj\-' a^'/\r,uj)Q^Kdr = T{uj)/\uj\ , (4-8) 

J r_ Jri 

while the second term becomes — vr/4, so this equation then gives 
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However, this limiting behavior can only exist if the buoyancy frequency A^^ is also small, 
since it provides a lower limit to uj as indicated in equation ( p-2| ). In addition, it is seen 

1 /2 

that Tj ~ ^> j {r/h)QQf), so this limiting behavior of be reached only for very large 

values of n. 



4.3. Neutrally Buoyant Disks 

The vertical eigenfunction equation ( P-28| ) and eigenvalue equation ( |4-2|) simplify 
greatly if the buoyancy frequency A^^ (cx dA/drj) vanishes. Equation ( |2-14| ) shows that 



this will occur if the specific entropy s is only a function of density, giving the barotropic 
equation of state p = p{p, s) = p{p). In practice, this can usually only be achieved by the 
entropy being constant (in the z direction), such as within convective zones of stars. 

The relation between pressure and density in the equilibrium (unperturbed) disk is 
determined by the equation of state, the radiative and turbulent advective transport of 
energy, and the viscous energy generation rate. The last two poorly known functions are 
critical in establishing the vertical gradient of the specific entropy s. However, as shown 
by Nowak & Wagoner (1995), with the maximum optical depths of accretion disks much 
less than within stars, advection can never dominate the transfer of energy. Therefore, the 
vertical gradient of the entropy cannot vanish. We also do not expect that s = s{p), so the 
buoyancy frequency should be nonzero, although it could be relatively small. 

Nevertheless, we shall now assume that the buoyancy frequency vanishes, in order to 
obtain relatively simple results that can be used as a guide to the nature of the eigenvalues 
\E'j and eigenfrequencies anj- In the next section, we relax this assumption to small 
buoyancy frequencies. 



When dA/drj = 0, equation ( P-14| ) gives p{r,z) = Fp{r)p^ . From vertical hydrostatic 



equilibrium [equation ( p-18| )] one then obtains 

P=Pc(r)(^l-|j , ^0 = ^ = ^- (4-10) 
(We note that the corresponding expression for rj^ in NWS lacks the factor of 2.) With 



Q = 0, equation ( |4-2|) then gives the vertical eigenvalue relation 

,2 



' 2E2(A;„,) ""^ Oi^!,' ^ ' 



where E{k) is the complete elliptic integral of the second kind. Employing the 
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local approximation, equations (^4-7|) and ( 4-11 ) allow the determination of \E'j and 
^nj = o'nj + fn^m for large eigenfrequencies. 

A graphical solution is presented in Figure ^, obtained in the following manner. 
Substituting equation ([4-7|) for o"„j into the expression for knjiojnj) in equation (|^ 



iTD gives the function k^j = fi{^j;n,Ni,N2), where the dimensionless parameters 



Ni{a) = Krn/^± ~ 1 and N2{cs{rm,0)/c,a) = VV^/{K,m^±am) ~ hm/Rm- Then equation 
( [4-ll| ) can be written in the form 

(j + Srr'^, = /2(vl>,-; n, T, N,, . (4-12) 

Therefore, in Figure ^ we have plotted [for a specific choice of a, F, and Cs{rm, 0)] each side 
of equation ( [4-12| ) versus \E'j, for various values of j (left side), n (right side), and m. From 
their intersections we obtain the eigenvalues \E'j. Then from equation ( [4-7|) we obtain the 
eigenfrequencies anj- 

In Tables |^, 0, and |^ we present (for three choices of a as well as m) the eigenvalues 
and corresponding eigenfrequencies for radial mode numbers n < 1, even vertical mode 
numbers J < 1, and odd vertical mode numbers j < 0. For comparison, we also indicate the 
corresponding values of the maximum eigenfrequency am- 

The expression ( [4-ll| ) provides the following bounds on the eigenvalue, since £'(0) = 7r/2 
and E{1) = 1: 

2(r - l)(j + < < Irr^T - l)(j + 6^)' • (4-13) 

E{k) is also monotonic between and 1, so the graphical solution shows that for any 
allowed value of j there is exactly one eigenvalue \Efj. For large enough values of j, it is close 
to its upper bound: ^ |vr2(r - l)(j + drf + 0(1). 

In the low-frequency limit (for m = 0), a'^j <^ Kq, equation ([4-ll| ) gives 
^ [(F - l)/2]n^{j + 5^)^ while equations(|3) and give 

J + 6r \ /F-l^l/^ 



^nj — 



poo 

j{o) ^ j{o) ^ / ^^^^^^ _ (4.14) 

J Ti 



4.4. Nearly Neutrally Buoyant Disks 

We now consider the case of buoyancy frequencies A^^ which are finite, but still much 
smaller than the characteristic frequency f2. Equation ( [4-lCI| ) then is modified to the form 

V = Pc{r){l - r/V%')'/^[l - fiv)] , (4-15) 
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where /(O) = and \f{ri)\ <^ 1. The vertical eigenvalue equation 
simple form if we choose f{ri) = \fQrf, with /q a constant: 



takes a particularly 



(r-l)(j + 5r 



2^2 (fc 



2 ^{knjj 



(4-16) 



with the complete elliptic integral of the first kind. Note that this result represents 

the first-order correction to that of the neutrally-buoyant disk for any slowly- varying (even) 
function f{r]). 



5. CONCLUSIONS 

The most observable g-modes will be those with the lowest values of the radial index 
n, since the radial compressions and rarefactions of the disk gas density and temperature 
will tend to produce the greatest modulations of the total luminosity . The luminosity 
modulations of modes with large values of the vertical index j should be less suppressed. 
These low n (high-frequency) modes are also those for which the local approximation is 
most valid. Therefore we can use equation ( p-16| ) to estimate (taking x ~ 1) that the 



effective radial width Ar = r_|_ — r_ of the eigenfunction is of order 

Ar ~ \/hRm oc \/cs(r„,0) . (5-1) 

Recall that Rm is an effective radius of curvature of the function — uj'^, defined in Section 
3.2. We find that for m = 0, Ro/ro varies from 0.408 for a = 0.0 {vq = 8.000) to 0.660 for 
a = 0.998 (ro = 2.449). For m = 2, i?2/r2 varies from 0.0444 for a = 0.0 (rs = 6.141) to 
0.00880 for a = 0.998 (ra = 1.242). 



We can also use equations (|2-25|) and (|2-26|) to obtain the following estimates 



of the relative magnitudes of the components of the fluid diplacement vector ^, in a 
local orthonormal basis, in terms of the scalar potentials 6V = K-K; = ^p/ (pP^) and 
W={k'^- uj^)-^dVr/dr: 

We also note that there are n radial nodes (corresponding to Unj) in W{r), j + 1 vertical 
nodes (corresponding to \l/j) in even Vr,{ri,r) for 11/9 < F < 7/5 and in odd V^(r7,r) for 
9/7 < F < 5/3, and j vertical nodes in even K)(r7, r) for 7/5 < F < 3 and in odd K)(^, for 
5/3 < F < oo. Only vertical nodes in the interval < rj < rjo have been counted. Since the 
pressure perturbation 6p oc K-K? oc {dW/dr)Vr^, it will instead have n + 1 radial nodes. 
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The accretion disk model chosen for Figures |^ and |^ and for Tables |l|, 0, and |^ is 
radiation-dominated in the region where the mode is trapped, corresponding to F = 4/3. 
There are thus j + 1 vertical nodes in both even and odd vertical eigenfunctions. From 
the results in these tables, we therefore see that the lowest even and odd temperature 
perturbations 5T oc 5p have one radial node. On the other hand, 5T has no vertical node 
for the lowest odd mode. The lowest even mode has one vertical node for low values of a 
and m, and no vertical node for the higher values. (Note, however, that the total actual 
number of vertical nodes will be twice as many for the even modes and one plus twice as 
many for the odd modes.) 

The eigenfrequencies — cr„j of the low n modes lie just below the maximum frequency 
(Tm- Using equations ( [4-7| ) and ( ^-111 ), the definition Rm = i^m/V^, and the approximation 
"ifjfl^/uj^j ^ 1 (which is seen to be reasonably valid from Tables the frequency 

separation is found to be approximated by 



C^m I '^nj I 



2F 



1/2 



>+l/2) 
(j + -^r) 



h{rr, 



n-Rr. 



(5-3) 



For the accretion disk model chosen for Figures ^ and ^ and Tables corresponding to a 
typical disk luminosity, we see from these tables that the eigenfrequencies \<Jnj\ of the lowest 
modes are indeed very close to the maximum g-mode frequency a^- This small splitting 
corresponds to a very thin disk [ h{rm)/rm = 6.6 x 10~^ for a = m = ], as indicated in 
equation (|5-3|) . 



Therefore the low-n g-mode inertial frequencies / = —a/2TT (now in usual units) of 
thin disks are essentially a function of only the mass and angular momentum of the black 
hole and the axial index m. For instance, we obtain 

/(m = 0) ^ (cVG'M)Ko/27r = {Q47TV2)-\cyGM)F{a) = 7U{MQ/M)F{a) Hz , (5-4) 

where F{0) = 1 and F{a) is the upper plot of Figure ^. It is seen that F{amax) — 3.443, 
if CLmax — 0.998. In order to determine both M and a, one needs to observe another 
eigenfrequency with a different dependence on a. 

For instance, in principle one could determine a via the frequency ratio 

/(m ^ 0)//(m = 0) = am/ (To = (k„ + mfi„)/Ko , (5-5) 

which is shown in the lower plot of Figure ^ for m = 1,2,3. Unfortunately, these ratios 
are close to their a = values unless a is fairly close to unity. And of course, any m ^ 
mode will not be observable if the disk is viewed face-on, and in general should produce less 
luminosity modulation (via Doppler boosting, etc.) than the m = mode. 
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Finally, we note that for a ^ 0.1 these results will also apply for accretion disks around 
slowly rotating neutron stars, since their external metric is Kerr through first order in a. 
However, the requirement that the disk be unaffected means that the neutron star must 
also be weakly magnetized {B <^ 10 gauss) and compact (soft equation of state). Fully 
relativitic numerical calculations ( Friedman & Ipser 1992; Cook, Shapiro, & Teukolsky 
1994) indicate that the maximum value of a = 0.6 — 0.8 over a wide range of equations of 
state. Typically, this maximum does not occur at the maximum (uniform) angular velocity. 

In future papers, we will extend this relativistic analysis in at least three directions: 

(a) We will investigate the characteristics of the other classes of modes, extending 
previous analyses of acoustic and corrugation modes. The dependence of the eigenfrequencies 
on the angular momentum parameter a is especially important, as indicated above. For 
example, the observation of three eigenfrequencies which have significantly different 
dependences on a would allow one to test the validity of the Kerr metric. 

(b) Numerical solutions for the eigenf unctions and eigenvalues will be obtained, in a 
manner similar to that employed in the pseudo-Newtonian analyses of NW1,2,3. Initially, 
the radial-WKB separated equations (|2-27|) and ( p-28| ) will be employed. 

(c) The assumption of thin disks will be relaxed. This will introduce many 
complications; such as noncircular flow, inclusion of more terms in the metric, and the 
probable need for a fully two-dimensional integration of the Ipser-Lindblom master 
perturbation equation. 

This research was supported in part by NASA grant NAG 5-3102 to R.V.W., NASA 
grant NAS8-39225 to Gravity Probe B, and a National Defense Science and Engineering 
Graduate Fellowship to D.E.L. We thank Lee Lindblom and Michael Nowak for useful 
discussions. 
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Table 1: Eigenvalues and eigenfrequencies of the lowest g-modes, calculated within the 
local approximation in neutrally buoyant disks. The angular momentum parameter a — 0, 
and the accretion disk model is the same as specified in Figure 5. For comparison, the 
maximum frequencies are (Tq = 0.0220971, ai = 0.0772040, and (72 = 0.1413801. As usual, 
the frequencies are expressed in units of c^/ GM. 



-24- 



Tfi 


71j 


9 

J 






— IT 


n 


n 


n 






n (^^9870^4 


n 

u 


1 


n 






n n'^9'^791 Q 


n 

u 


n 


1 

X 






n n'^9QQ81 Q 

U-UOZjC/C/OXC/ 


n 


1 


1 




1 .\JOOOO X 


n (^^975971 


n 


n 


_1 


nrlrl 


U.U X(JUO(J 


n n'^94fi4'iS 


n 

yj 


1 

± 


_1 

X 


nrlrl 


U.UO / c/TiX 


03193175 











odd 


4.419339 


0.03295735 





1 





odd 


4.429199 


0.03263055 


1 

± 


n 

u 


n 

u 


t; Veil 


9 409449 


n 1914(S9qi 


1 

± 


1 

± 


n 




9 4n'^41 7 


n 191*^8^97 

U. XZi XOOtJZi / 


1 

X 


n 


1 

± 




R 1 9^97'i 


U. XZi X'4:OUOX 


1 

± 


1 

X 


1 

X 


c veil 




19143896 


1 

± 


n 


_1 

~ X 


nc\c\ 


n 78fi847 


n 1914*^909 


1 

1 


1 


1 

1 




U. ( oozou 


D 1919Q97Q 




n 


n 
u 




4 )^^n9nfi 




J- 


1 

J. 







'-t.oOyjijO X 


n 19149098 


2 





-1 


even 


0.065171 


0.22445146 


2 


1 


-1 


even 


0.065678 


0.22432569 


2 








even 


2.499892 


0.22450645 


2 


1 





even 


2.500037 


0.22448958 


2 





1 


even 


8.244807 


0.22451026 


2 


1 


1 


even 


8.244903 


0.22450099 


2 





-1 


odd 


0.865544 


0.22450043 


2 


1 


-1 


odd 


0.865750 


0.22447153 


2 








odd 


4.960238 


0.22450891 


2 


1 





odd 


4.960353 


0.22449696 



Table 2: Same as previous table, for angular momentum parameter a — 0.500. For 
comparison, the maximum frequencies are Uq — 0.0331210, (7i = 0.1215017, and (72 = 
0.2245149. 
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Table 3: Same as previous table, for angular momentum parameter a — 0.998. For 
comparison, the maximum frequencies are ctq = 0.0760771, ui — 0.4277195, and (T2 — 
0.8456764. 
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Fig. 1. — The radial dependence of the square of the fundamental free-particle frequencies 
that govern the modes of the disk: Keplerian and radial (k) and vertical {fi±) epicyclic. 
Three values of the black hole angular momentum parameter a — cJ/GM^ are chosen. 




Fig. 2. — The functions which determine the radii r-{a) and r_|_((T) between which 
the eigenfunctions are concentrated. Also indicated is the maximum value of the 
eigenfrequency \a\, and the inner radius n of the disk. For this plot, we choose a — 0.5 
and m — 1 (with the case m — corresponding to use of the dashed horizontal lines) . 




Fig. 3. — The radii between which a g-mode eigenfunction is concentrated is plotted versus 
its eigenfrequency. 




Fig. 4. — (a) The dependence of the maximum radial (m = 0) eigenfrequency on the black 
hole angular momentum parameter a = cJ/GM'^, relative to its value at a = 0. (b) The 
ratio of the maximum eigenfrequency of some higher m modes to that of the radial mode. 
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Fig. 5. — (a) The black hole angular momentum dependence of the radius to which 
r_ and r+ converge as \a\ — > cr^. Also shown is the radius of the inner edge of the 
disk, (b) The dependence of the fractional effective width of the lowest eigenfunction, 
^f^/f^m = [^"+(0") ~ ''^-{'^)]/f^m, on the angular momentum of the black hole. The same 
values of m are chosen as in (a). The accretion disk model is specified by F = 4/3, a locally 
isentropic equation of state, and speed of sound corresponding to a luminosity L = O.lLEdd 
from a radiation-dominated disk. 




Fig. 6. — Graphical determination of the lowest eigenvalues assuming the local 
approximation for a neutrally buoyant disk. The accretion disk model is the same one 
used in the previous figure, and we have chosen a = 0.5. The curves are the right-hand-side 
of equation (4-12) for n — 0,1,2 (difference invisible) and the values of m indicated, while 
the straight lines are the left-hand-side of equation (4-12). 



